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Abstract. An indefinite stochastic Riccati Equation is a matrix- valued, highly non- 
linear backward stochastic differential equation together with an algebraic, matrix 
positive definiteness constraint. We introduce a new approach to solve a class of such 
equations (including the existence of solutions) driven by one-dimensional Brownian 



motion. The idea is to replace the original equation by a system of BSDEs (with- 
er 

out involving any algebraic constraint) whose existence of solutions automatically 

enforces the original algebraic constraint to be satisfied. 
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1 Introduction 

Stochastic matrix Riccati equations were first introduced by Bismut [1] in his study of some 
stochastic control problems. A very special class of these equations is the so-called quadratic 
backward stochastic differential equation (BSDE). The existence and uniqueness of solutions 
for such BSDEs remain a largely open problem, particularly for BSDE systems; but see [6], [5], 
[2] and the references therein for recent progress. The indefnite stochastic Riccati equations 
(SRE) were first formulated in [3], motivated by the introduction of the indefinite stochastic 
linear-quadratic control problems. Such an equation is typically matrix- valued, highly nonlinear 
(not even quadratic), and involves a positive semidefinite constraint in addition to the backward 
equation. In pE], the uniqueness of solutions to the SRE was established in the greatest generality 
based on a control argument, but the existence was solved only for several very special cases. 
The general existence remains to this date a significant open problem. 
The SRE is a BSDE over a running time interval [0,T]: 



dP = J2 A J dWJ 
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subject to the constraint that 



K = R + J2 D 'j PD J > ° 

3=1 



XD 
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in the matrix sense, and subject to the terminal condition that P{T) = H which is Tt- 
measurable. In this formulation, the time parameter t is omitted for simplicity, the capital 
letters A, B,C, D,A, Q and P are real matrix valued (adapted) processes and D' means the 
transpose of D etc. All the matrix processes are square with the same dimension n, P and A,, 
j = 1, 2, • • • , k, are unknowns and all the other parameters are given, and W = (W l , • • • , W k ) 
is k- dimensional standard Brownian motion. 



The given matrix R in defining K is called the gauge matrix, which is an adapted process. 
The SRE fll.ip - fll.2D is indefinite, if the gauge matrix R is allowed to be indefinite, i.e., R can 
have zero or negative eigenvalues. 

The problem is to look for square integrable adapted processes P and (Aj) satisfying the 
corresponding stochastic integral equations as well as the constraint (11.21) . Moreover, in view 
of the proved uniqueness of solutions, a solution matrix P must be symmetric as long as the 
parameters Q, R and H are symmetric. 

We believe that the existence for the general SRE (ll.ip - (ll.2l) in high dimensions will remain 
to be an open question for some time. The main challenge, apart from the highly nonlinear 
nature of the BSDE (II. ip and the fact that the equation is matrix-valued, stems from the 
presence of an additional algebraic constraint (11.21) . In this note, we develop a new approach to 
solve a class of SREs driven by one-dimensional Brownian motion. The main idea is to consider 
a system of BSDEs satisfied by (K, K~ x ) without any algebraic constraint. It turns out that 
the existence of solutions to the equation satisfied by K~ x can be established independently, 
which in turn will ensure the validity of the original constraint K > 0. 

2 The main result 

In this paper we consider the SRE driven by a one- dimensional Brownian motion, with the 
matrix process D being invertible. 

Therefore, in the remainder of this note, W is a one- dimensional standard Brownian motion 
on a complete probability space (Q, J 7 , P), and (J-i)t>o is the Brownian filtration generated by 
W. Without lose of generality, we may assume that D = I, and we study the following BSDE 

dP = AdW - [PA + MP + C'PC + AC + C'A + Q] dt, 

+ [PB + C'P + A}K- 1 [B'P + PC + A}dt, te[0,T}, (2.1) 

P(T)= H, 



where T > is fixed throughout the paper and the random symmetric matrix H is bounded 
and J-r-measurable, subject to the constraint that 

K = R + P>0. (2.2) 

For simplicity we assume that the coefficients A, B, C, Q are bounded, square n x n matrix 
valued adapted processes, and in addition Q is symmetric. 

We are interested in the indefinite gauge case, that is, R may have zero or negative eigenval- 
ues, not necessarily being positive definite. The most interesting case in applications is when 
R is a matrix-valued Ito process; so we assume that it has the representation 

R(t) = R(0) + f Fds+ f GdW (2.3) 

Jo Jo 

where F, G are bounded, adapted and measurable symmetric matrix- valued processes. 

By a solution (P, A) we mean a pair of (J-^-adapted, measurable, square integrable and 
matrix-valued processes P = (P(t))te[o,T\ and A = (A(t))te[o,T\ such that K(t) = R(t) + P(t) > 
in the matrix sense for all t G [0, T], a.s., and 



P(t) = H- J AdW 

+ f [PA + A'P + C'PC + AC + C'A + Q] ds 
- J [PB + C'P + A] R- 1 [B'P + PC + A] ds 



(2.4) 



for t G [0, T], a.s., where the stochastic integral is understood in the Ito sense. A solution (P, A) 
is called bounded if P is bounded. Clearly, if (P, A) is a solution, then P must be a continuous 
matrix-valued semimartingale. 

Due to the presence of K~ l in (12. ip . it is natural to rewrite (12.11) in terms of K = R + P 



and A = A + G. This can be achieved by making substitutions in (12. ip : P by K — R and A by 
A — G, leading to the following BSDE 



dK = AdW- 



KA + A'K + Q 



dt + 



KB+A-R 



K- 



B'K + K-R! 



K(T) = R{T) + H, 



dt, 



(2.5) 



where 



Q = Q + F + C'RC + R (BC -A) + (C'B' - A') R 



(2.6) 



and 



A = A-BC, R = RB + C'R + G. 



We are now in a position to state our main result. 



Theorem 2.1 Assume 



(2.7) 



R = RB + C'R + G = 



(2.8) 



and 



Q = Q + F + C'RC + R (BC -A) + (C'B' -A')R> 0, 



(2.9) 



(ii) G and F are bounded adapted measurable processes such that R(T) + H > 0, and there is 
a constant 5 > such that (R{T) + H)~ > SF 

Then there is a unique solution (P, A) to the SRE A2.1\) - fi£Jfy . Moreover, P + R is bounded. 

The uniqueness has been established in j3], Theorem 3.2. The existence, on the other hand, 
is known in the so-called definite case, namely, when R > 0, Q > 0, H > 0; see [7] (where there 
is an additional assumption that D = 0). The existence when all these matrices are allowed to 
be indefinite is investigated in [1] for several special cases. It should be noted that the existence 



of (12. lp is by no means unconditional; the problem is to find sufficient conditions under which 
the existence holds. One of the conditions, (12. 81) . of Theorem 12. II stipulates that R satisfies the 
following Ito equation 

dR = -(RB + C'R)dW + Fdt (2.10) 

where F can be arbitrary (up to the required Lebesgue integrability) . The other condition, 
(j2.9p . requires an "overall" positive semidefiniteness in place of that of individual matrices. 
We will make further comments on the conditions of the preceeding theorem in Section [5j 



3 A linear equation 



We need a result about the representation for a linear matrix-valued BSDE. Consider 



dY = J2 U i dWJ 

3=1 

Y T = H , 



m K 

YA + A'Y + J2 Hy^ + ^ {u j C j + Cfy + Cj-YC,-) + Q 

»=1 .7 = 1 



dt, 

(3.1) 



where A, C, E and Q are n x n matrix valued, adapted and bounded, H is bounded and 
Tt- measurable, and W is a standard Brownian motion of dimension k. The BSDE is linear so 
there is a unique solution. Choose a standard Brownian motion W of dimension m, which is 
independent of W. For any p 6 M™ and < t < T, let £ be the solution to the linear stochastic 
differential equation 



d$ = a^s + J2 c£ dW3 + J2 Eitdtir, Ct = p- 



3=1 



»=1 



(3.2) 



Lemma 3.1 Under the above assumptions and notations, we have 



p'Y(t)p = E 



&H&+ / CQZsds 



F, 



(3.3) 



Proof. Applying Ito's formula to 



J(s) = e s Y(s)C s = (Y(s)C s ,Cs) 



(3.4) 



we obtain 



k m 

d,J = -£'Q£ds + Z'J2 fa + Cj Y + YCj) dW J C + £' J2 ( Y ^ + ^ Y ) d ^- ( 3 ' 5 ) 

j=l i=l 

Integrating from t to T, and conditional on Tt we obtain ( I3.3p . ■ 

Lemma 3.2 // m addition H > 51 for some constant 5 > and (J > 0, then the solution Y 
to A3. 1\) satisfies 

Y(t)>8e-M r - t) I VtG[0,T], a.s. } (3.6) 



where 



(3 = ess sup {-2 inf f'i(s,u;)f,0}. 

wen,s<T ICI =1 

Proof. Let p6R" and £ solve (13. 2D . Applying Ito's formula to |£| 2 to obtain 



(3.7) 



d|£| 2 = 2eA^ds + J2\C^\ 2 ds + J2\E^\ 2 ds 



i=i 



*=1 



+2 J] gCjtdW* + 2J2 C'E^dW 1 . 

3=1 



t=l 



Integrating over [£, T] and taking conditional expectation on J 1 ^ we obtain 

/s I" / fe m \ 

> b| 2 -/?y E[|e r | 2 |J- 4 ]rfr. 



^ 



The Gronwall inequality then yields 

E[|£r| 2 |^]>b|V^-<). 
Finally, it follows from Lemma [3. II that 

p'Y(t)p > E l&Hfr JF t l > 5\p\ 2 e-^ T - l) 
which implies (13. 6p . ■ 



4 Proof of Theorem 12.1 



The remainder of the paper is devoted to the proof of Theorem 12.11 

To handle the positive defmiteness constraint (I2.2p . we couple the BSDE (12 .5p together with 
another BSDE for X = K~ x , and consider the resulting system of BSDEs without the explicit 
constraint K > 0. This last constraint will be implied by the existence of solutions to this 
system of BSDEs. 

Therefore we next derive the BSDE for X = K~ l which can be obtained from the identities 
XK = KX = I. In fact, by integrating by parts, 

dX = -X(dK)X - (dX, dK)X. (4.1) 

In particular the martingale part of X is — XAXdW; so (dX, dK) = —XAXAdt. Substituting 
this equation into (14. ip to obtain 

dX = -X(dK)X + XAXAXdt . (4.2) 



Using (12. 5p we obtain a BSDE that X = K l must satisfy, that is 



dX 



-XAXdW + X 



KA + A'K + Q 



Xdt 



-X 



KB-R+A 



X 



B'K -R' + A 



Xdt + XAXAXdt. 



(4.3) 



Setting Z = —XAX, and using the fact that KX = XK = I, we obtain its equation 



dX 



ZdW 



AX + XA' - BXB' + BZ + ZB' 



dt 



ZR'X + XRZ 



dt 



+ 



XQX + BXR'X + XRXB' - XRXR'X 



dt, 



(4.4) 



X(T)= (R(T) + H)-\ 



Notice that all the terms involving A have been canceled out thanks to the assumption that 
the driving noise W is one dimensional; so (14.41) no longer contains A. This reveals another 
significant feature of the SRE, that is, equation (14.41) for the inverse matrix X = K~ l is 
itself closed, in the sense that it does not depend on K or A. Therefore we can solve (14.41) 
independently without the prior knowledge that X is the inverse of K. Let us call (14.41) the 
inverse equation associated with the SRE (I2.ip - (12.2I) . 

Therefore, if we are able to solve (14.41) to get (X, Z) with X > on [0, T], then (P, A), where 
P = X- x -R and A = -X^ZX-G, is a solution to (O). In particular, R+P = K = X~ x > 
is satisfied automatically. 

Now we return to the BSDE (12.51) for K and we wish to rewrite it in terms of (Z, X). There 
are several ways to do this because of the relations XK = KX = I, and we will choose one 
which will serve our propose in this paper. In (12.51) replace K~ l by X and replace XA by —ZK 
to obtain 



dK = AdW 



KA + A'K + Q 



dt 



KB+A-R 



XR'dt 



+ 



KB+A-R 



{XB 1 - Z) Kdt, (4.5) 

K(T) = R(T) + H. 

We consider (14.41) and (14.51) together as a single system of BSDEs, and ignore the fact that 
X is the inverse matrix of K as well as the constraint K > 0. This system can be solved one by 
one: we can solve (14.41) first to obtain (X, Z), and then solve (14. 5p regarding (X, Z) as known 
parameters. This is actually the approach we will follow. 

Under our assumption that R = 0, our basic BSDEs (\AA\\ and ( 14.51) are significantly simpli- 
fied. In fact 



dX = ZdW 



AX + XA' - BXB' + BZ + ZB' 



dt + XQXdt, X(T) = (R(T) + H)-\ (4.6) 



and 



dK = AdW - 



KA + A'K + Q dt + (KB + A) (XB' - Z) Kdt, K(T) = R(T) + H. (4.7) 



Lemma 4.1 Suppose ((X,K), (Z, A)) is a bounded solution to l[4-bV and j^.7| ). Then XK 
KX = I. 



Proof. Let Y = KX — I. Then applying Ito's formula to ( 14. 6 p and ( 14.71) we obtain 



dY = K(dX) + (dK)X + {dK, dX) 

= UdW + UB'dt + YA'dt - A'Ydt + YQXdt 

+(AXB' - KBZ - AZ)Ydt 
= UdW + UB'dt + Y (A' + Qx)dt 

+ 



AXB' -A' -(KB + A)Z 



Ydt, 



10 



where U = KZ+AX. This is a linear BSDE with the terminal value Y(T) = K{T)X{T)-I = 0. 
The uniqueness of solution to the linear BSDE then yields Y = 0. ■ 

BSDE ( 14. 6 p is matrix valued with a quadratic term in the drift. If Q > 0, then it is a special 
case of a definite SRE whose solvability has been established by Bismut [T] and Peng [7J . In our 
case where Q > 0, we use an approximation scheme, adapted from [TJ, to prove the existence 

of TO). 

Lemma 4.2 Let n be a bounded, n x n symmetric matrix-valued J-'x-'m^asurable random vari- 
able. Then there is a unique adapted bounded solution (X, Z) to the BSDE 

dX = ZdW + (AX + XA' - BXB' + BZ + ZB') dt + XQXdt, 
X(T) = rj. (4.8) 

If in addition n > 0, then X(t) > for all t G [0, T\. 

Proof. To show the existence of the BSDE (I4.8p . we consider the following iteration 

dX (n+i) = Z {n+i) dw + \AX {n+1) + X (n+1) A' - BX {n+1) B' + BZ (n+1) + Z {n+1) B'] dt 
+X {n+l) QX {n) dt + X {n) QX {n+1) dt - X {n) QX (n) dt, 
X {n+1) (T) = 7] 

which is a linear BSDE, whose unique solution defines (X^ n+1 \ Z^ n+1 ^). Since r/ > 0, each 
X^ > (Lemma EH). 

Let Y&) = iW - X(" +1 ) and U^ = Z^ - Z^ n+1 \ Then the pair (Y< n ), U^) satisfies the 
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following stochastic equation: 

dY in) = jj(n) dw + \AY {n) + Y {n) A' - BY {n) B' + BU in) + U {n) B' 
+Y {n) QX {n) dt + X^QY^dt 



d,t 



+ 

Y {n) {T) = 0. 



X (n) QX (n - 1} + X {n - l) QX {n) - X {n) QX {n) - X {n - 1] QX {n - 



i) 



d* 



Note that, since Q > 0, the symmetric matrix 

xWqj(»-i) + x(«-i)qjW_j(«)qiW_X(«-i)qi(«-i) = _( I (")_ I («-i))q( I («)_ J ("-i)) < 

for each n, hence Y( n > > (Lemma 13.11) . It follows that X*™) > and X( n ) is decreasing in 
matrix sense, and therefore has a unique limit denoted by X. It is then routine to show that 
{ZW : n > 1} converges as well (in L 2 ([0, T] x ft, dt ® dP)) to a limit process Z. Then (X, Z) 
solves flUHD - 



Since Q > and 77 > in matrix sense, X(t) > for all t G [0, T]. Furthermore, if Q > 
and ?7 > 0, then X(t) > for any < t < T. To see this, we apply the representation ( 13. 31) to 
A; = 1, m = 0, A = A + ^XQ and B = C\. Fix t <T and define £ s by solving the corresponding 
SDE (13.21) . If p ^ 0, then, by the uniqueness of linear BSDE for (13. 2 p with terminal £t, we can 
conclude that £r 7^ 0. As a result, £' t H£t > 0, and therefore p'X(t)p > a.s.. ■ 

Lemma 4.3 Under the same assumption as in Lemma \4^ and if in addition rj > SI for some 
5 > 0, then X" 1 is bounded. 

Proof. Since X is bounded, we conclude that 

A, = ess sup {- inf (2A£, + XQ£, £), 0} 

is finite. It then follows from Lemma [3.21 that X(t) > 5e _/3oT J so that X -1 is bounded. ■ 
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Lemma 4.4 Let n = {R{T) + H) > 51 in Lemma \4-S\ and (X, Z) be the corresponding 



solution. Then P = X' 1 - R solves the SRE fO)-/ TO)) . 

Proof. Let K = X~ l which have been proved to be a bounded matrix valued semimartin- 
gale. Let A = —KZK. Since KX = XK = I, it follows from Lemma [4. II that (K,A) satisfies 
(14.51) . Therefore (P, A), where P = X^ 1 — R and A = A — G, in turn solves (12.11) . Moreover, 
K = X^ 1 > 0, namely, the constraint (j2.2p is satisfied. ■ 

The uniqueness for SRE is well known. The proof of Theorem 12.11 is complete. 

5 Discussions and Examples 

In this section we discuss about the assumptions of Theorem 12.11 and give examples for illus- 
tration. 

First of all, if R = 0, then condition (12.81) holds automatically, and ( 12. 9 j) is equivalent to 
Q > 0. In addition, assumption (ii) in Theorem 12.11 boils down to H > 0. (The condition 
H -1 > 51 is implied by H > and the fact that H is bounded.) In this case, our result 
improves Theorem 5.2 in j3], since here we do not need to assume C = and H~ l is bounded. 

If R is a constant yet indefinite matrix, and B — C — 0, then again ( 12.81) is satisfied, and 
(12. 9p reduces to Q — RA — A'R > 0. In this case we recover Theorem 5.3 of |4J. However, 
from Theorem 12.11 we immediately realize that the assumption B = C = is far from being 
essential. Indeed, in the case when R, B and C are non-random matrices, the essential condition 
is RB + C'R = 0, which can be satisfied easily by infinitely many non-zero matrices B and C 
and indefinite matrices R. In this case, the condition Q — RA — A'R > should be replaced by 
Q + R + C'RC + R (BC - A) + (C'B' - A')R> 0, where R denotes the derivative of R, which 
is zero if R is a constant matrix. 

As a matter of fact, we can "generate" many generally indefinite, adapted processes R 
satisfying condition (12.81) . To see this, let S be the solution of the following matrix- valued, 
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sample-wise ODE 



dS 



C'W t?„bw ni cd mil 



Fe BW _ C , SB _ (p^ s _ SB < 



dt 



with any given initial state, where F is any given adapted process so that e c Fe BW is inte- 
grable over t e [0, T] a.s.. Define R = e~ c w Se~ BW . Then Ito's formula yields that R satisfies 
( |2~TU]) . namely, (J23J holds. 

If condition (12. 8 ft does not hold (i.e. R ^ 0), then, we will need to study the general inverse 
equation, (14. 4p . This is a very interesting BSDE, since it is matrix- valued involving a cubic 
term of X. In general, its global existence is not guaranteed. For example, suppose that R= I, 
A = 0, B = 0, and that the terminal and Q are non-random. Then (j4.4p becomes 



dX 



XQX - XXX 



dt 



whose solution may explode in a finite time (and thus the corresponding SRE can not have 
a global solution). An interesting and challenging open problem is to identify the "weakest" 
condition on R so that the cubic BSDE (|4.4p admits a solution. 

On the other hand, the condition (R(T) + H) > 61 for some 5 > is more technical 
than essential. One can weaken this condition by incorporating more involved technicalities 
in our analysis. However, the main goal of this note is to introduce and highlight the main 
approach, that is to use a system of BSDEs to substitute the original mix of a BSDE and 
an algebraic constraint, to solving the indefinite SRE. Therefore, we have preferred not to let 
undue technicalities distract the main idea. 

Finally, let us remark that the proof of Lemma I4.2[ along with Lemma I4.4[ has indeed 
suggested a numerical scheme to solve the indefinite SRE. 
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